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ABSTRACT 

We report results from a deep Jansky Very Large Array (JVLA) search for 12 CO 
J = 1 — line emission from galaxies in a candidate galaxy cluster at z ~ 1.55 
in the COSMOS field. We target four galaxies with optical spectroscopic redshifts 
in the range z = 1.47 — 1.59, consistent with the likely redshift for the candidate 
galaxy cluster. Two of these four galaxies, ID 51613 and ID 51813, are nominally 
detected in CO 1 — line emission at the 3 — 4cr level. We find CO luminosities of 
(2.42 ±0.58) x 10 10 K km s" 1 pc 2 and (1.26 ±0.38) x 10 10 K km s" 1 pc 2 , respectively. 
Taking advantage from the clustering and expanded 2-GHz bandwidth of the JVLA, 
we perform a search for emission lines in the proximity of optical sources within the 
field of view of our observations (60"). We limit our search to galaxies with A"s < 23.5 
(AB) and z phot = 1.2 - 1.8. We find 2 bright optical galaxies, ID 51207 and ID 51380, 
to be associated with significant emission line peaks (> 4cr) in the data cube, which 
we identify with the CO 1 — line emission. To test the reliability of the line peaks 
found, we performed a parallel search for line peaks using a Bayesian inference method. 
Both CO line emitting candidates are identified with probabilities of 13% and 72% 
that there are line emitting sources in each case, respectively. Monte Carlo simulations 
show that such associations are statistically significant, with probabilities of chance 
association of 3.5% and 10.7% for ID 51207 and ID 51380, respectively. Modeling of 
their optical/IR spectral energy distributions (SED) indicates that the CO detected 
galaxies and candidates have stellar masses and star formation rates (SFRs) in the 
range (0.3 — 1.1) x 10 11 Mq and 60 — 160 Mq yr" 1 , with star formation efficiencies 
(SFEs) comparable to that found in other star-forming galaxies at similar redshifts. 
By comparing the space density of CO emitters derived from our observations with 
the space density derived from previous CO detections at z ~ 1.5, and with semi- 
analytic predictions for the CO luminosity function, we suggest that the latter tend to 
underestimate the number of CO galaxies detected at high-redshift. Finally, we argue 
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1 INTRODUCTION 

One of the major goals of modern observational cosmology is 
to understand how the gas in the diffuse interstellar medium 
(ISM) of galaxies converts into stars and how both phases 
evolve with cosmic time. A major advance has been the 
determination of the evolution of the SFR density of the 
Universe. Deep optical and radio surveys indicated that the 
SFR density steadily increases from z = 5, with a peak at 
z ~ 3 — 1, and steeply declines from z = 1 to the present 



(e.g., Lillv et al.l 19961; Madau et al.ll996l 
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Steidel et all 19991; Giavalisco et al.lbo04l; 


Ouchi et al.ll2004 


Bouwens et al.l2007l; Smolcic et al.l2009l). Although the con- 



tribution from luminous, merger-driven starburst galaxies 
at high-redshift to the SFR density appears to be signifi- 
cant, the population that drives such evolution appears to 
be formed by more quiescent massive star -forming galaxies 
l|Daddi et alj|2007l ; iRodighiero et al.l 120111 ). The close rela- 
tionship between the SFR and the amount of molecular gas 
in galaxies, from which stars form, suggest that the evolu- 
tion of the SFR density is the result of the evolution of the 
molecular gas density in galaxies across cosmic times. 

Measuring the evolution of this molecular gas den- 
sity is thus critical to understand the formation of molec- 
ular gas and stars in galaxies. For this, it is neces- 
sary to perform blind surveys of the sky that allow 
us to measure the molecular gas content in a CO flux 
limited sample as a function of redshift and luminos- 
ity (e.g. observations in the local Universe and simula- 
tions by Keres et al.l 20031 ; lObreschkow fc Rawlings 2009; 
lObreschkow et all 12009 ) . The technical limitations (e.g. 
bandwidth, sensitivity) of current submillimeter and radio 
facilities, along with the intrinsic faintness of the CO emis- 
sion lines used to measure the amount of molecular gas in 
galaxies have, however, precluded such studies. The advent 
of sensitive int erferometers such a s the Jansky Very Large 
Array (JVLA; iPerlev et al. 2011T) and the Atacama L arge 
Millimeter Array (ALMA; IWootten fc Thompsonll2009r i will 
make these studies possible. 

A huge step toward this goal has been the recent de- 
tection of significant CO line e mission in massive star- 
forming galaxies at high redshift dDaddi et al1l2008l . [2010a; 
ITacconi et al"1l2010l ; ICeach et all 1201 ll ). The large molecu- 
lar gas masses derived are comparable to that observed in 
submillimeter galaxies (SMGs), in the range ~ 10 10-11 Mq, 
however these galaxies are forming stars at ~ 5 — 10 times 
lower rates, with typical SFRs in the range ~ 50 — 300 Mg 
yr _1 . These findings imply that they have low SFEs and 
thus consume the gas at longer timescale s, in agreement wit h 
their long duty cycle times of > 0.5 Gyr (|Daddi et al.ll2007l ). 
From the analysis of the CO molecular gas and Ha ionized 
gas kinematics, it is expected that between 1/3 to 2/3 of the 
star fo rming galaxies at z ~ 2 are c o nsistent with rotatin g 
disks i|Forster Schreiber et all 120091 ; ITacconi et all l2010h . 
and suggest that the dominant mechanism to cool gas into 
molecular clouds is d riven by cold gas accretion rather tha n 
galaxy interactions (|Keres et all 120051 ; iDekel et all 12009^ . 
Follow-up multi-transition CO observations recently discov- 
ered that a cold gas component, with gas excitation condi- 
tions similar to th at seen in local spiral gal a xies, is present 
in th ese objects (|Dannerbauer et al.l 120091 ; lAravena et al.l 
20 !<: :■. Remarkably, observations at z — 0, z = 0.5 and 



z — \ — 2 suggest that there has been a strong evolu- 
tion in the molecular gas fractions in star forming galaxies 
over the past 10 G yr (|Tacconi et al.ll2010l ; ILagos et al.ll201ll ; 
iGeach et al1l201lh . 

In this paper, we present results of a deep search for 
CO molecular gas from galaxies that are part of a candidate 
galaxy cluster at z ~ 1.5 with the JVLA. The scope of this 
paper is to perform the first attempt to conduct an efficient 
blind search for CO emission from galaxies located in a part 
of the sky that likely contains a higher than average number 
of CO emitters. The clustering of galaxies in angular scales 
and along the line of sight makes it possible to effectively in- 
crease the number of galaxies observed within a single point- 
ing and frequency tuning with the JVLA. We use the ex- 
tensive multi-wavelength dataset available for the COSMOS 
field, including the accurate photometric redshifts, to inform 
the location of optical galaxies and to guide the search for 
CO emission. In section 2, we present the properties of the 
(overdense) target field and the details of the JVLA obser- 
vations. In section 3, we present the main results from this 
paper. We report the detection of CO 1 — emission in two 
spectroscopically confirmed galaxies and analyze the 45.5 
GHz continuum emission and properties of a spectroscopi- 
cally confirmed radio loud galaxy. In sections 3.3 and 3.4, 
we present the search for line emission peaks from optical 
galaxies in the field. Significant line peaks associated with 
optical sources that have compatible photometric redshifts 
are identified as CO line candidates. The line search is per- 
formed using two different algorithms. From sections 3.5 to 
3.7, we analyze the likelihood that such associations are real 
based on both a statistical approach and on the detection 
of the CO candidates in the IR wavelength regime. In sec- 
tion 4, we discuss the derived SFEs for the two CO detected 
spectroscopically confirmed galaxies and for the two CO can- 
didate emitters, as well as the implications of these results 
for future experiments to measure the CO source counts at 
high redshift. We adopt a concordance ACDM cosmology 
throughout with Hp = 71 km s" 1 Mpc" 1 , Q M = 0.27, and 
n A = 0.73 (|Spergel et al.ll2007l ). 



2 OBSERVATIONS 
2.1 The COSMOS field 

COSMOS is the largest deep survey carried out with the 
Hubble Space Telesco pe (HST) covering ~ 2 deg 2 in the sky 
|Scoville et all 120071 ). Extensive imaging of the COSMOS 
field has been performed from the X-rays to the radio wave- 
lengths. This includes complete optical/near-IR coverage in 
broad and intermediate bands with several ground-based 
observatories including the Subaru telescope, the Canada 
France Hawaii Telescope (CFHT), the United Kingdom In- 
frared Telescope (UKIRT), and the Kitt-Peak National Ob- 
servatory (KPNO); IR imaging with the Infrared Array 
Camera (IRAC) and the Multiband Imaging Photometer 
for Spitzer (MIPS) on board of the Spitzer Space Telescope; 
and radio imaging with the VLA at 20 cm. Deep imaging 
at far-IR wavelengths with the Herschel Space Observatory 
has been obtained as part of the Photodetector Array Cam- 
era and Spectrometer (PACS) Evolutionary Probe (PEP) 
and the Hersche l Multi-tiered Extragalactic Survey (Her- 
MES) programs (|Lutz et all 1201 ll ; lOliver et all 120101 ). This 
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Figure 1. (Left:) Bi'K color composite of a 3' X 3' region around the candidate cluster center. Cyan contours show the XMM Newton 
X-ray emission at 2,3 and 4c significance. White circles show spectroscopically confirmed galaxies with redshifts 2 S pec = 1.470 — 1.595. 
Yellow circles represent spectroscopically confirmed galaxies with z Bpcc = 1.2 — 1.8. A large 30" radius circle indicates the location and 
field of view of our deep JVLA observations of the molecular gas. (Right:) Close-up to the central 80" X 80" around the JVLA pointing 
position. Labels indicate the sources ID tags. 



Table 1. Observed CO properties for galaxies with spectroscopic redshift 
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a J2000 b 
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10 02 40.47 


+01 34 41.2 


1.523 




< 0.053 


< 0.7 




51130+ 


10 02 42.25 


+01 34 32.5 


1.519 




< 0.044 


< 0.6 




51613 


10 02 43.36 


+01 34 20.9 


1.516 


1.517 


0.20 + 0.05 


2.42 + 0.58 


8.7 + 2.1 


51858 


10 02 40.43 


+01 34 13.1 


1.560 


1.556 


0.10 + 0.03 


1.26 + 0.38 


4.5 + 1.4 



t Radio galaxy; a COSMOS ID; 6 Position of the optical source; c Optical redshift; d CO redshift; e Spatially and velocity integrated 
line flux; ' CO luminosity; 9 Gas mass obtained using chqo = 3.6 Mq (K km s -1 pc 2 ) -1 



has been complemented by an intensive spectroscopic cam- 
paign with gr ound-based observatori es as part of the zCOS- 
MOS survey (|Lillv et all 120071 . 120091 ). From the rich multi- 
wavelength imaging, accurate photometric redshifts were 
computed using 31 UV/optical/IR bands (|Salvato et alJ 
120091 ; lllbert et all 120091 : (Salvato et ai1l201ll ). These reach 
uncertainties of Sz ~ 0.2 — 0.3 at z ~ 1.5 (see below). For 
details in the optical/IR im a ging and c atalogs of the COS - 
MOS field, see lCapak et~ai1 (j2007h and lllbert et~afl l|2009l ). 
The 1.4 GHz radio observations in this region of the COS- 
MOS field reach an rms level of 30 jj,3y beam - , and hence 
only sources with fluxes > 100 ^iJy beam -1 can be reli- 
ably detec ted. A complete descrip t ion of the radio im aging 
is given in ISchinnerer et al.l (|2007l ). iBondi et all (|2008l ) and 
ISchinnerer et all I 20ld ). The targeted region lies right out- 
side the region of the COSMOS field defined by the HST 
imaging coverage, and therefore it is not covered by Ad- 
vanced Camera for Surveys (ACS) F814W band imaging. 
Hence studies of the optical morphologies of galaxies is not 
possible. 



2.2 Our target: Discovery of a candidate galaxy 
cluster at z ~ 1.5 with the red-sequence 
technique 

In this section, we present the identification of a substantial 
candidate galaxy cluster at z ~ 1.55 based on four pieces of 
evidence: (i) the detection of a significant (7x) overdensity 
of galaxies with photometric redshifts z — 1.5 — 1.6; (ii) the 
tentative detection of X-ray emission; (iii) the formation of 
a red sequence of a few passive galaxies in the field; and (iv) 
the existence of a radio loud galaxy with a confirmed spec- 
troscopic redshift consistent with the redshift of the overden- 
sity. We remark that regardless of whether this constitutes 
a virialized structure (e.g. a galaxy cluster) or not, the over- 
density of galaxies within a confined redshift range allowed 
us to observe the CO emission from 4 spectroscopically con- 
firmed massive galaxies in that redshift range (Fig. [TJ, and 
to search for blind CO line detections from optical galaxies 
in the field. 

Accurate photometric redshifts computed with 31 
bands indicate that 26 galaxies with K < 23 mag (AB) 
and within ~ 3' from the central X-ray position (see below), 
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Figure 2. Projected number density of galaxies with z p hot = 
1.5 — 1.6 and SFR > 5 Mq yr _1 in a 3' X 3' region around the 
cluster center. The grayscale and contours represent the density 
of galaxies in this field, &c, given in terms of the average density 
So of similarly selected galaxies in the COSMOS field, 5 = Sq/Sq. 
Contours range from 2 to 7 in steps of +1. The red circle repre- 
sents the location of the JVLA pointing and primary beam (PB) 
FWHM. 
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Figure 4. The distribution of photometric redshifts of galaxies 
located within r2oo from the center of the cluster candidate. The 
open histogram shows all the galaxies in the field. The red his- 
togram shows the redshift of the red-sequence identified galaxies, 
and the dashed histogram indicate the redshifts of galaxies shown 
in blue in Fig. |3] 




Figure 3. Color-magnitude diagram (z — [3.6] vs. [3.6]) of galax- 
ies located within r2oo (0.9') from the center of the cluster can- 
didate. The gray dashed line shows a mo del red sequen c e for a 
cluster formation redshift of z f =3 from Li dman et alj J200ST> . 
The gray star symbol represents the the characteristic magnitude 
m* for passively evolving galaxies. Blue filled squares correspond 
to galaxies in the redshift range z = 1.5 — 1.7. The red filled circles 
represent galaxies located in the same redshift range but identi- 
fied to belong to the forming red-sequence. Small black circles 
show background/foreground galaxies in the field. 



have 2 p hot = 1-5 — 1.6. Based on their best-fit spectral en- 
ergy distributions (SED), we find that 16 of the 26 galaxies 
are consistent with significant star-forming activity and SFR 

> 5 Mq yr _1 . For galaxies in this redshift range, this repre- 
sents an overdensity of ~ 7x compared to field star-forming 
galaxies at this level. This is possibly the most compelling 
evidence of the existence of a cluster in this region. To illus- 
trate this, we constructed a number density map of galaxies 
with photometric redshifts in the range 1.5 — 1.6 and SFR 

> 5 Mq yr _1 in this region, using a procedure similar to 
the one introduced by iDresslerl (|l980l ). For this, we created 
a grid of 300 x 300 positions spaced by 2" centered at the 
inferred cluster X-ray position, and computed the density 
from the distance to the lOth-nearest neighbour galaxy, dio, 
to each grid point. The number density is thereby computed 
as raio = 10/(7rdi ). Similar approaches have be en applied 
to rich galaxy clusters (e.g., iGuzzo et al.l 120071 '). Figure [2] 
shows the density map obtained, convolved with a Gaussian 
kernel that has a full-width at half maximum (FWHM) of 
9". The number density, S, is given in terms of the aver- 
age density of field galaxies with z p hot = 1.5 — 1.6 and SFR 

> 5 Mq yr" 1 in the COSMOS field, S = 0.7 arcmin" 2 . 
The source density in this region is ~ 7x larger than this 
average. In terms of the standard deviation in the distribu- 
tion of source density in the field, ao = 0.24 arcmin -2 , this 
corresponds to a fluctuation of ~ 18ct with respect to the 
average source density in the field. The average and stan- 
dard deviation values were obtained by creating a density 
map of galaxies with z p hot = 1-5 — 1.6 and SFR > 5 Mq 
yr _1 in a nearby ~ 20' x 20' region that did not include the 
overdensity, using the same procedure as above. From the 
distribution of density values in this map, the average and 
standard deviation were computed. 

The structure is tentatively detected in the XMM New- 
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ton X-ray images at the 3<r level (Fig. [TJ, and a bright, 
extended radio source close to the center of the cluster 
is seen in the VLA 1.4 GHz images (see section 3.2). 
This radio source was spectroscopically confirmed to be at 
2 = 1.5185 ± 0.002, based on the detection of the Mgll, 
[NeV], [Oil] and [Nelll] emission lines. The ubiquity of 
radio sources in high-reds hift dense large-scale structures 
l|Milev fc De Breuck 120081 ). lends support to the existence 
of this galaxy cluster candidate. Details of this source are 
given in section 3.2. 

The cluster red-sequence technique represents an effec- 
tive method to select galaxy clusters. It is based on the fact 
that galaxy clusters typically contain a uniform population 
of early-type galaxies. When using two photometric band 
filters to bracket the location of the 4000 A break, which 
is typically a strong feature in passively evolving galaxies, 
early-type galaxies can be identified as the brightest and 
reddest objects in the field, forming a tight relation in color- 
magnitude diagrams. Such a "red-sequence" is pivoted by 
the brightest cluster galaxy, typically has a small scatter, 
and has been shown to hold out to z ~ 1.5 (|Lidman et alJ 
|200S| ; iBielbv et al.ll20ldl ; iTanaka et al.ll20ld ). 

Based on the red-sequence te c hnique modifi e d for 
high-redshift dLidman et al.l 120081 : IWilson et aD 120091 : 
iFinoeuenov et aT 2010h . we found evidence for the forma- 
tion of a red sequence of passive galaxies in this field (Fig. 
[l]). We used the (z — [3.6]) galaxy colors, which select clus- 
ters at z > 1.2. Figure [3] shows the color-magnitude diagram 
obtained for galaxies located within a radius of 0.9' of the 
cluster center (0.5 Mpc at z = 1.55). When selecting objects 
with photometric redshifts in the range z p hot ~ 1.5 — 1.7, 
a number of bright, red galaxies compared to the bulk of 
"blue" galaxies is visible. Figure [4] shows the distribution 
of photometric redshifts for galaxies in the field, highlight- 
ing the distribution formed from the red and blue galaxies 
from Fig. [3] From the redshift distribution for all galaxies 
in this region (empty histogram), it can be seen that there 
is a redshift peak in the distribution at z ~ 1.5 — 1.6. 

In addition to the available optical spectroscopic red- 
shifts provided by zCOSMOS, a dedicated spectroscopy 
campaign targeting galaxies within a 3' x 3' region around 
the center of the cluster candidate was performed using the 
Deep Imaging and Multi-Object Spectrograph (DEIMOS) 
on the Keck telescope and Multi-Object Infrared Camera 
and Spectrograph (MOIRCS) on the Subaru telescope (Fig. 
[1] Tanaka et al. in preparation). Due to bad weather condi- 
tions, however, it was not possible to obtain redshifts of the 
passive galaxies in the field. Nine sources within this region 
were confirmed to have spectroscopic redshifts in the range 
1.47 — 1.59, which corresponds approximately to the range 
covered by the JVLA observations (see next section; Fig. Q]), 
including the galaxy associated with bright 1.4 GHz radio 
emission. This redshift range represents a large range in ve- 
locity, > 10000 km s _1 , which could point to a filamentary 
structure nature for the overdensity rather than a collapsed 
one. 

Based on the spectroscopic and photometric redshift 
evidence, we pointed our JVLA observations targeting four 
spectroscopically confirmed galaxies that were located close 
to the maximum galaxy density in this field. In order to 
observe the CO 1—0 emission from the four spectroscopically 
confirmed galaxies, which includes the bright radio source, 



we tuned the JVLA receivers to cover the redshift range 
z ~ 1.47 — 1.59 (see next section). Hereafter, and in the 
absence of spectroscopic evidence confirming this structure, 
we refer to this as a cluster candidate or overdensity. We 
stress that the existence of an overdensity in photometric 
redshift and angular scales, added to the large bandwidth of 
the JVLA, enables us to search for CO 1 — emission from 
galaxies in this field irrespective of whether this overdensity 
is a cluster or not. 



2.3 JVLA observations 

Observations of the 12 CO J = 1 — emission line (u rcB t = 
115.271 GHz) in the cluster galaxies at z ~ 1.5 were taken as 
part of the JVLA Residence Shared Risk Observing (RSRO) 
program (program ID: AA330). The target field was ob- 
served during March 2010 to January 2011 and centered 
at ajaooo = 10 h 02 m 41.6 s , <5 J2 ooo = +01°34'24."o. The ob- 
servations were made using the Q-band receivers (tunable 
range: 40-50 GHz) and taking advantage of the expanded 
bandwidth of the WIDAR correlator (» 2 GHz). 

Based on early spectroscopic and photometric redshift 
estimates that indicated a cluster redshift of z w 1.53, we 
tuned our observing frequency to 45.500 GHz. We used a 
wide bandwidth mode, with 2 intermediate frequency (IF) 
bands (or basebands: AC and BD), each formed by 8 spec- 
tral windows (SPW). Each SPWs was composed of 256 chan- 
nels with a resolution of 2 MHz per channel. The SPW IDs 
range from to 15. Combined, they cover a total bandwidth 
of 2048 MHz (~ 2 GHz). Since each SPW loses sensitivity 
at the edge channels and because at the time when the ob- 
servations were taken the WIDAR correlator could not be 
configured to overlap different SPWs, we used two 2048 MHz 
bandwidth spectral setups shifted by 5 channels from each 
other. Two channels at each edge of the SPWs were thus 
flagged. This avoids spectral gaps, obtaining an even fre- 
quency coverage. At the edge SPW of each IF, the passband 
deteriorated more than in the other SPWs. Thus, 5, 3, 4 and 
4 additional channels at the external IF borders of SPW ID 
0, 7, 8, and 15, respectively, were flagged. This results in 
a coarser frequency coverage and a spectral gap of 8 chan- 
nels (16 MHz or 105 km s _1 ) in the middle of the band. 
Such problem was unforeseen prior or during the observa- 
tions. Hence, the observations covered the frequency ranges 
44.486 - 45.492 GHz (AC) and 45.508 - 46.514 GHz (BD), 
or approximately 1 GHz per band. This implies that we are 
able to detect the CO 1—0 line emission from galaxies in 
the redshift range z = 1.478 — 1.591 with a negligible gap of 
Az « 0.0009. 

Although for detection experiments the D-array con- 
figuration is better suited, the observations were taken in 
C-array configuration due to scheduling constraints. This 
implies that sources extended by more than 1" were slightly 
resolved with the consequent loss in sensitivity. The nearby 
source J1038+0512 (~ 0.5 - 0.7 Jy; ~ 9° distance to tar- 
get) was observed every ~ 10 min for gain calibration. The 
brighter sources J1256-0547 (3C279; ~ 16-21 Jy) and 
J1331+3030 (3C286; « 1.4 Jy) were used for bandpass and 
flux calibration, respectively. Pointing calibration was per- 
formed on J1038+0512 every ~ 55 min using the C-band 
receivers at 4.9 GHz. 

Data calibration was done using the Common Astron- 
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Figure 5. JVLA CO 1 — emission line intensity map (PB corrected), averaged over 230 and 500 km s — 1 and centered at 45.810 
GHz and 45.107 GHz for ID 51613 and ID 51858, respectively. The background corresponds to a Bi + K color-composite image. The CO 
emission is shown in contours in steps of ler, starting at ±2cr, with <r = 120 and 40 /^Jy beam , for ID 51613 and ID 51858, respectively. 
The corresponding beam sizes are 2.1" X 1.9" and 1.1" X 1.0", respectively. Dashed contours represent negative signal. In the case of ID 
51613 the resolution has been downgraded to increase the signal to noise. 



0.6 



0.4 



0.2 



0.0 



-0.2 



■ ■ ■ i ■ ■ ■ ■ i ■ ■ ■ ■ i ■ ■ 1 ■ i ■ ■ ■ ■ i 
ID 51613 



z DPT =1.516 



lu 




\ n41 



-1000 -500 500 1000 
Velocity (km/s) 



^ 0.2 




-1000 -500 500 1000 
Velocity (km/s) 



Figure 6. JVLA 45.5 GHz close-up spectrum towards the optical galaxies ID 51613 and ID 51858 at 93 and 132 km s per channel 
resolution, respectively. The emission profiles were measured at the peak position of the integrated emission close to the optical source. 
The arrow represents the frequency that corresponds to the CO 1 — line at the optical spectroscopic redshift. The emission features 
identified with the CO 1-0 line emission are highlighted by a Gaussian fit (blue solid line). 



omy Software Application (CASA) package. Time ranges 
with poor visibilities as well as the edge channels where 
the bandpass deteriorated, as explained above, were flagged. 
The final integration time amounts to ~ 20 hrs on-source. 
At the observing frequency, the JVLA primary beam (PB) 
FWHM is ~ 60". To obtain flux measurements, we decon- 
volved the visibilities using the CLEAN task with natural 
weighting and a further taper of 0.8", and applied the cor- 
responding primary beam correction. We thus obtained a 
final resolution of 1.1" x 1.0" and an average rms of 60 /iJy 
beam - per 100 km s -1 channel. This corresponds to a phys- 
ical scale of ~ 8.5 kpc at the redshift of the overdensity and 
thus well suited to extract their extended emission without 
losing many long baselines. 



3 RESULTS 

3.1 CO 1 in two spectroscopically confirmed 
galaxies 

Our JVLA observations were designed to cover the CO 
1 — line emission from four galaxies that have an ac- 
curate optical/near-IR spectroscopic redshift in the range 
z = 1.478 — 1.591 and located spatially within the JVLA 
PB FWHM. 

Since the frequency of the CO 1 — line emission could 
be slightly offset from the frequency implied by the optical 
redshift, we first searched for a CO emission peak in the 
spectra along the line of sight. From this, we estimated the 
velocity integrated line flux emission at the CO peak fre- 
quency, if any. If no emission was detected in the spectra 
at the source position, we measured a limit to the emission 
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Figure 7. Radio continuum emission towards the massive radio 
loud galaxy ID 51130. A Bi + K$ color composite image is shown 
in the background. The VLA 1.4 GHz emission is shown in yellow 
contours with levels in steps of 2cr, starting at ±2cr, with <r = 
30 /^Jy beam -1 . The 45.5 GHz continuum emission is represented 
by white contours in steps of la = 10 ^tjy beam -1 , starting at 
±2(T. Dashed contours represent negative fluctuations. 



based on the optical redshift assuming a typical line FWHM 
of 300 km s -1 . 

Two out of the 4 galaxies with optical spectroscopic 
redshift, ID 51858 and ID 51613, were nominally detected 
in CO 1 — emission with significances of ~ 3<r and ~ 4<r, 
respectively. Figures [5] and [6] show the velocity averaged CO 
emission maps and the CO line profiles, respectively, com- 
puted at the position of the CO image peak emission. Table[T] 
summarizes the positions, redshifts, physical properties and 
derived parameters for these sources. 

After inspection of the CO images with 1.1" x 1.0" res- 
olution, the emission from ID 51613 appeared to be slightly 
resolved along the major axis. In order to perform flux mea- 
surements in this case, we applied a taper of 1.7" to the 
visibilities, leading to an image resolution of 2.15" x 1.90". 

From a Gaussian fit to the spectral line profiles, we mea- 
sured formal line widths of 360 ± 220 km s" 1 and 200 ± 80 
km s -1 at central frequencies of 45.1087 GHz and 45.8030 
GHz for ID 51858 and ID 51613, respectively. Spatially inte- 
grated flux densities were obtained using a two-dimensional 
Gaussian fit to the CO maps. These results imply CO red- 
shifts of 1.5554±0.0008 and 1.5167±0.0004 for ID 51858 and 
ID 51613, respectively. This leads to a difference of ~ 590 
km s _1 with respect to the optical redshift in the former 
case. Such difference is large compared to that observed in 
some o f the massive disk galaxies reported bv lDaddi et all 
l^OlOal ), of ~ 200 - 300 km s" 1 , however it is within the un- 
certainties implied by the low significance of this detection 
and the accuracy of the optical spectroscopic redshift deter- 
mination, typically Sz ~ 0.001, and well within the typical 
velocity differences between CO observations and rest-frame 
UV for SMGs and QS Os, with a typical range - 500 km 
s _1 l|Greve et al.ll2005h . 



Figure 8. Optical to mid-IR SED of ID 51130. Gray filled circles 
with error bars represent the multi-wavelength photometry and 
the solid line shows the best fit template at the source's spectro- 
scopic redshift z = 1.519. 
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Figure 9. Comparison between spectroscopic and photometric 
redshifts for galaxies in the COSMOS field at z = 1.5 — 1.6. A 
dispersion a = 0.275 is measured from the mean /i = —0.008. 



3.2 A massive radio galaxy at z = 1.519 

Figure shows the VLA 1.4 GHz image of a prominent 
double peaked radio source identified in the field (see also 
Fig.QJ. The radio source is associated with the bright optical 
galaxy ID 51130, which has an optical spectroscopic redshift 
of 1.519 (Table[T]). Both 1.4 GHz peaks are identified as the 
source's radio lobes. 

Based on the presence of radio lobes, a strong contribu- 
tion from AGN activity would be expected. Since the source 
is not detected in the XMM Newton X-ray maps, we ex- 
plored the shape of the mid-IR SED since AGN sources tend 
to have a steeply rising SED from 3.6 /im to 24 /im. For this, 
we compare the mid-IR color s of ID 511 3 with three differ- 
ent co lor-co lor criteria fr om Lacy et ail (|2004l ). IStern et al.l 
(|2005h and llvison et all (|2004T i. ID 51130 has flux densi- 
ties S3.6 M m = 87 fljy, S4.5(im = 97 jLtJy, Ss.Sum = 84 /xJy, 

Sg.o M m = 78 /iJy and SW^m = 220 /^Jy, respectively, be- 
ing clear it has a relatively flat mid-IR SED (Fig. [HJ. 
While, the criterion proposed by Lacy et al. selects AGN 
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based roughly on the colors log(5 , 5.8nm/5 , 3.6 M m)> —0.1 and 
log(5 , 8.0nm/S , 4.5 M m)> —0.2, Stern et al. uses limits roughly 
defined by 0.7 < (ms.8 — ms.o) < 1.6 and (7713.6—7714.5) > 0.3. 
Similarly, Ivison et al. defines a region occupied by AGN 
with Ss.opm/Si.s/im > 1 and SW/zm/S'g.opm < 5- Based on 
this, we find that ID 51130 only satisfies the Lacy et al. 
AGN criterion. The source is located outside the AGN re- 
gions defined by Stern et al. and Ivison et al., however in 
both cases this could still be consistent with a Seyfert 2 
galaxy. From this, we conclude that there is not a substan- 
tial contribution from an AGN to the galaxy's emission at 
optical/IR wavelengths. 

Based on the lack of evidence pointing to an important 
AGN contribution, we fitted a pure galaxy model (with no 
AGN contribution) to the source's photometry (Fig. [5] see 
section 3.7). Note that fitting pure galaxy templates to ob- 
jects with no substantial AGN contrib ution such as S e yfert 2 
galaxies is consistent with results from lSalvato et all i 201ll ). 
The host galaxy's SED is best fitted by an 0.5 Gyr old galaxy 
template with a SFR of 25 Mq yr~ that appears to have 
already formed a good fraction of its stars, as suggested by 
the large stellar mass (Table 3). This is consistent with the 
non-detection of CO line emission from this object down to 
a limit of L' co = 0.6 x 10 10 K km s" 1 pc 2 (3a). Since the 
nature of the host galaxy is unknown, we assume a CO lu- 
minosity to gas mass conversion factor aco = 3.6 Mq (K 
km s -1 pc 2 ) -1 , which is the typical value f ound for disk 
galax ies at high-redshift (|Daddi et alj|2010al ; iMagdis et al.l 
l201ll ). and derive a limit to the gas mass of this object of 
2 x 10 10 Mq (3a). We notice that this value of aco traces 
the total amount of gas, including H2, Helium and HI, with 
the vast majority of the hy drogen (> 90%) in molecular 
form (see Dad di et al]|2010al ). Using a conversion factor of 
~ 1 Mq (K km s _1 pc 2 ) -1 would only put a stringent but 
less conservative limit. This limit is also consistent with the 
molecular gas masses typically f ound for radio AGN in th e 
nearby Universe (e.g. Table 4 in ISmolcic fc Riechersll201lT ). 

To measure the continuum emission of this object at 
45.5 GHz, we averaged the whole 2 GHz data set over chan- 
nels and created a continuum image, inverting the visibilities 
using natural weighting with an additional tapering of 0.8". 
Figure [7] compares the optical view of this galaxy with the 
radio emission observed at 1.4 GHz and at 45.5 GHz. The 1.4 
GHz radio maps reach an rms level of ~ 30 /xJy beam -1 in 
this region of the COSMOS field. Two positive blobs, close 
to where the bulk of 1.4 GHz emission is located, can be 
distinguished in the 45.5 GHz continuum images, although 
at very low significance, ~ 2 — 3a. By fitting a double com- 
ponent two-dimensional Gaussian to the 45.5 GHz map, we 
find an integrated flux density of 145 ± 50 pJy and 23 ± 8 
/xjy for the western and eastern radio peaks, respectively. At 
1.4 GHz, both peaks are detected with in tegrated fluxes of 
400±30 /-tJy and 520±9 /jJy, respectively ijSchinnerer et all 
l2007l ; iBondi et al.ll2008l ; ISchinnerer et alj|2010h . Deeper ob- 
servations at 45.5 GHz would be needed to properly measure 
the spectral index of the radio lobes. 



3.3 Search for CO line emission in galaxies 

We took advantage of the unprecedented multi-wavelength 
coverage of the COSMOS field, as well as the accuracy 
of photometric redshifts computed with 31 UV/optical/IR 



bands filbert et all 120091 ; ISalvato et all |201ll ). to identify 
possible CO emission line peaks from optical galaxies within 
the field of view (60"). 

We selected galaxies with Ks < 23.5 mag (AB) and 
photometric redshifts in the range 1.2— 1.8. The magnitude 
limit is set to guarantee that the photometric redshift de- 
termination is accurate, based on the comparison between 
spectroscopic and photometric redshifts for the COSMOS 
field; and that we trace galaxies with stellar masses > 2 x 10 9 
Mq. As mentioned above, our JVLA spectroscopic observa- 
tions trace the CO 1 — emission line in the redshift range 
« 1.48 — 1.59. In this range, the photometric redshift accu- 
racy is Az = |zs pcc — -Zphotl ~ 0.275, with a small system- 
atic offset Zspcc ~ Zphot +0.008 (Fig. [9)|. Hence, to comprise 
most of the galaxies in the redshift range 1.48 — 1.59, we 
conservatively selected galaxies with photometric redshifts 
in the range 1.2 — 1.8. Based on the spectroscopic redshifts, 
we expect 5-10% of galaxies to be photometric redshift out- 
liers. From this, we have that a total of 26 galaxies are thus 
selected, of which up to 3 are expected to lie outside the 
specified redshift range. 

For each selected optical source, we extracted a 10" x 
10" cutout from the full JVLA data cube, and searched 
for significant (> 4a) emission line peaks lying within a 
1.0" radius from the optical source position. We performed 
the search making use of the AIPS task SERCH. For this, 
we used a data cube that had been de-convolved using a 
Briggs weighting scheme, with a robust parameter of 1.0, 
and an additional tapering of 0.8", without PB correction. 
The weighting scheme without PB correction ensures a sim- 
ilar noise behavior across the image, for all source positions. 
The SERCH task uses a Gaussian kernel to convolve the 
data cube along the frequency axis with an expected input 
line width, and reports all channels and pixels having a sig- 
nal to noise ratio over the specified limit. We experimented 
with various Gaussian kernel line widths, from ~ 200 to 500 
km s -1 and found that an optimal width is using ~ 350 
km s . Using larger line widths does not increase the num- 
ber of coincidences between optical sources and line peaks, 
however smaller widths (~ 200 km s -1 ) tend to produce a 
significantly larger number of line peaks (30% more) with 
no optical associations, which are likely spurious. 

A similar approach to search for CO candidates at high 
redshift is to blindly search for peaks in the full field-of- 
view of the JVLA observations, and then cross-correlate 
the possible detections with the optical catalog, using the 
same constraints on magnitudes and photometric redshifts 
and searching for nearby associations. Since such approach 
is still limited by the existence of optical sources, it essen- 
tially leads to the same results. We focused only on line 
peaks associated to relatively bright optical sources since 
this approach will produce more reliable results. Note that 
the frequency coverage of our observations also enables us 
to detect CO 2 — 1 emission from galaxies in the redshift 
range z = 3.956 — 4.182, however such sources are expected 
to be very faint in the optical, being more difficult to find a 
clear optical counterpart identification. Search for line peaks 
associated to z ~ 4 sources is out of the scope of this paper. 

We find that 2 optical sources, ID 51207 and ID 51380, 
can be identified with emission line peaks from the data 
cube. The spatial coincidence in both cases is better than 
0.6". Table [2] summarizes the results and Figs 1101 and 1111 
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Table 2. Observed line properties for CO-identified galaxies with photometric redshift in the range 1.2 — 1.8 
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51207 10 ft 02 m 43.400 s +01°34'33.12" 0.60 1.529 1.530 0.085 ± 0.021 1.03 ±0.25 3.7 

51380 10 h 02 m 40.810 s +01°34'27.51" 0.34 1.450 1.551 0.090 ± 0.021 1.12 ±0.26 4.0 



a COSMOS ID; 6 Position of the optical source; c Optical photometric redshift; d CO redshift; e Spatially and velocity integrated line 
flux; ' CO luminosity; 9 Gas mass obtained using a galactic conversion factor c*co = 3-6 Mq (K km s -1 pc 2 ) - 1 . 



show the extracted spectra and images. In both cases, the 
line peaks have significances > 4a, however in the case of ID 
51207, the positive line peak in the spectrum is followed by 
a region of noise channels with negative flux density. This is 
likely because this source is located at the border of the con- 
sidered JVLA PB FWHM coverage, where the noise starts 
to increase. 

An additional line association is found for an optical 
source for which a secure optical spectroscopic redshift of 
1.240 (ID 51861) had been obtained. Its optical redshift dis- 
cards the possibility of it being a CO detection, since CO 
1 — could only be detected in the range « 1.48 — 1.59. 
At the redshift of this source, no bright emission line from 
molecules is expected. Therefore, we regard this as a spuri- 
ous detection. This identification is discussed in more detail 
in the next sections. 

None of the CO sources with optical spectroscopic red- 
shift were found using this method. This is because, at the 
resolution of the searched data cube (~ 1"), ID 51613 has a 
peak significance of only ~ 3.2a. In this case, it was neces- 
sary to deconvolve the images to a lower resolution of ~ 2" 
to maximize the signal-to-noise. In the case of ID 51858, the 
significance of the CO detection is too low to be picked up 
by our signal-to-noise limit. 

3.4 Line search from Bayesian inference 

We performed an independent search for line peaks us- 
ing model fitting to the data cube upon the princi- 
ples of Bayesian infere nce with the MultiNEST algorithm 
l|Feroz fcH obson 200j|). This method provides an efficient 
way of sampling from the posterior probability distribution 
and calculating evidence values. We use a galaxy model char- 
acterized by a set of parameters O (e.g. x,y positions, line 
profile central frequency), such that spatially a source is de- 
scribed by a circular or elliptical Gaussian, and in the fre- 
quency domain the rotation curve is approximated using 
either a Gaus sian line profile or Brandt's parametrization 
l|Branddl 19601 1. A full de scription of this approa ch to source 
line searches is given bv lFeroz fc Hobsonl l|2008l ). 

Selection between two models Ho (that there is no 
source present) and Hi (that there is a source present) 
can be performed by comparing their evidence values, Zi 
and Zi, via R = ZijZ%. The probability that there is a 
source present in the data is then computed from the rela- 
tion P = R/(l + R). 

The four galaxies that have optical spectroscopic red- 
shifts in the range covered by the JVLA observations were 
searched within a radius of 1" of the optical identification, 
and central frequencies accurate to within 1 channel width 



(~100 km s -1 ) from their optical spectroscopic redshifts. No 
line detection is found for ID 51130 and ID 51480, however 
we find 55% and 41% probability that there is line emission 
from ID 51613 and ID 51858, respectively. Probability val- 
ues > 50% are considered as good evidence that there are 
line detections, given that at this level < 30% of the iden- 
tified peaks are expected to be spurious line detections. In 
the case of ID 51858, the lower probability is due to the sig- 
nificant offset between the measured line frequency and the 
frequency implied by the optical spectroscopic redshift. 

The remaining galaxies with photometric redshifts were 
searched within a radius of 1" from the optical identifica- 
tions, but central frequencies covering the full bandwidth of 
the JVLA observations. Using this method, the two identi- 
fied CO source candidates presented in Section 3.3, ID 51207 
and ID 51380, were found to have probabilities of 13% and 
70% that there is a source present, respectively. ID 51861, 
however, has a probability of 53%, even though it is a con- 
firmed spurious detection. Such values must therefore be 
taken with caution, as they give an estimate of the reliabil- 
ity that the line peak can be described with a physical model 
of how a CO emission line typically looks like, and they do 
not take into account the properties of the identified optical 
sources. The low probability measured for ID 51207 is likely 
due to the fact that this source lies at a channel, an position 
(border of the image), where the noise degrades. 

Finally, we remark that a total of 7 line candidate de- 
tections with P > 50% are found using this method, 6 of 
which are also found with SERCH. Three of them corre- 
spond to sources studied in this work: ID 51380, ID 51613, 
ID 51861. The rest, despite their high probability, are peaks 
with lower signal to noise values, as measured directly in the 
JVLA maps, and therefore not considered further. 

3.5 Probability of chance association 

We estimated the probability that a line peak would be as- 
sociated by chance within 1.0" from an optical galaxy by 
performing a Monte Carlo simulation of the distribution of 
significant line peaks in the JVLA field. 

Before starting with simulations, we first carried out 
a blind search for all positive peaks that have a signal to 
noise ratio > 4.0 within the JVLA data cube, using a Gaus- 
sian kernel width of 350 km s _1 . This search yielded a total 
of 61 line peak candidates within the 60" JVLA primary 
beam. The peaks appear to be slightly clustered in some re- 
gions of the image and thus they are not strictly random. 
From the 61 peaks, 3 correspond to the optical identifica- 
tions presented in Section 3.3, of which one is spurious as 
noted in the previous section, and 9 constitute additional 
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Figure 10. JVLA CO 1—0 emission line intensity map (PB corrected), averaged over 480 and 260 for ID 51207 and ID 51380, respectively. 
The background image corresponds to a Bi + K color-composite. The candidate CO emission is shown in contours in steps of lcr, starting 
at ±2cr, with a = 53 and 57 /xjy beam -1 , respectively. 
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Figure 11. JVLA 45.5 GHz spectrum towards ID 51207 and ID 51380. Channel separation is 16 MHz or ~ 105 km s per channel. 
The emission profiles were measured at the peak position of the integrated emission. The identified emission feature is shown by a blue 
horizontal bar with a width equal to the line FWHM. 



line-peak/optical-source associations that do not meet all 
the specified Kab < 23.5 and z p hot = 1.2 — 1.8 criteria. All 
these extra associations coincide with optical sources in the 
redshift range 1.3 — 2.4. A total of 49 of the 61 line peaks 
(80%) do not have a clear nearby optical counterpart and 
likely correspond to spurious detections. This means that up 
to 20% of the cases could correspond to real CO line emit- 
ting galaxies. We note that we are biasing our counterpart 
selection by using only the optical i + -band selected source 
catalog. Thus, line peak detections with no optical counter- 
part could still represent real sources. For instance, highly 
obscured sources or at high-redshift (z > 4) are typically 
missed by optical catalogs and they can only be revealed in 
IRAC images (or at far-I R/submillimeter wavelengths; e.g., 
iDannerbauer et al.| [2008). As mentioned previously, search- 
ing for these kind of sources is out of the scope of this paper. 

To compute the probability P of chance association be- 
tween a random peak detection and an optical source within 
1" radius, we performed Monte Carlo simulations of line 
peaks in the JVLA region. In each of the 10 4 samplings, we 



generated 61 peaks drawn by the real spatial distribution of 
peaks in the field (we measured this by constructing a den- 
sity map of peaks using a 7th- nearest neighbor method), and 
thereby computed P as the fraction of simulated line peaks 
that are associated to our two optical sources: ID 51207, ID 
51380; and also to the spurious identification of ID 51861 
for comparison. 

Figure [12] shows the computed probability of chance as- 
sociation for these three sources (filled diagram), as derived 
from our simulations, compared to the probability that any 
of the other optical sources that meet our selection crite- 
ria would have been associated to a line peak by chance 
(solid histogram). For several optical sources that meet our 
selection criteria, including our line identified targets, the 
probability to find a line peak within 1" is non-negligible, 
but in all cases less than 20%. While for ID 51207 the prob- 
ability is small (~ 3%), in the cases of ID 51380 and ID 
51861 there is ~ 10% probability that they are associated 
by chance. Such probabilities represent an upper limit, since 
the actual matches in those cases were better than 1". This 
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Figure 12. Distribution of the probability P that an optical 
source is associated by chance within 1" with any of the 61 sig- 
nificant line peak in the observed region. The dashed histogram 
shows the value of P for the optical sources found to be associated 
with emission line peaks. The solid empty histogram shows P for 
the other optical sources, which met the selection criteria, if they 
were associated with significant line peaks. Given the low number 
of significant emission line peaks in the cube, an association by 
chance is unlikely. 



procedure, however, assumes that the optical selection crite- 
ria provides good galaxy candidates with redshifts accurate 
to within Sz < 0.2 in a case by case basis. 

In the case of the identification of the optical source 
ID 51861 with a line peak, it is important to note that 
even if the photometric redshift estimate (z p hot = 1.436) 
would have agreed better with the spectroscopic one (z spcc — 
1.240), we would have still selected this source as a possible 
line candidate given our selection range z p hot = 1.2 — 1.8. 
However, a critical difference between source ID 51861 with 
respect to ID 51207 and ID 51380 is that the two later cor- 
respond to faint IR detections (see section 4.1). Clearly, de- 
tection of IR emission, which is directly linked to dust and 
gas content, makes a CO detection more likely. The only 
way to avoid spurious detections is to impose tighter selec- 
tion criteria, including tighter photometric redshift ranges 
and detection of sources in the IR bands, at the expense of 
missing a larger number of emission line candidates. 

Finally, we estimated the false match probability by 
computing the number of associations between negative 
peaks in the data cube with SNR < —4a, searched in the 
same way than the positive peaks as above, with optical 
sources that meet our specified selection criteria. From a to- 
tal of 70 negative peaks, we find only one association within 
1" from an optical source, and thus imply a ~ 2% false 
match probability for our positive detections. 



3.6 Dust properties 

An important ingredient to characterize the properties of 
galaxies and to test the reliability of the CO candidate 
sources is their far-IR emission produced by dust. For the 
derived CO 1 — luminosities, it would be difficult to not 



detect the sources in the Spitzer and Herschel SPIRE maps 
for typical gas to dust ratios and SFRs. Detection in these 
IR bands and consistency of the expected IR luminosities 
can give strong supporting evidence of the reality of the CO 
line emission. 

A number of the galaxies targeted by the JVLA ob- 
servations have far-IR measurements provided by the Her- 
schel PACS and SPIRE imaging surveys of the COSMOS 
field (|Lutz et al.ll201ll ; lOliver et all 12012 ). Since the galax- 
ies in the targeted region are strongly clustered, the emis- 
sion is blended in the Herschel images. Fluxes from the 
PACS/SPIRE maps were thus extracted using the 24 / tm po - 
sitions as priors, as explained in detail bv lLutz et al.l (|201lh . 
This provides reliable measurements of the far-IR fluxes, for 
galaxies detected in the Spitzer MIPS 24 fim images. 

Within the JVLA field of view (Fig.[TJ , we find 10 MIPS 
24 /im sources, 8 of which have photometric redshifts > 1.0. 
Five of these 8 correspond to sources detected with SPIRE 
at 250 /im with z p hot = 1-4 — 1.6, 3 of which are also de- 
tected at 350 £tm with signal to noise ratio > 4a. Four of 
the 5 SPIRE detected galaxies have associated CO detec- 
tions. The remaining 3 out of the 8 high-redshift 24 /im 
sources, are either not detected by Herschel or lie at z > 2. 
None of the 8 high-redshift 24 ^m sources are detected with 
PACS at 100 or 160 /mi, or with SPIRE at 500 /im. These 
data, however, provide important constraints to the far-IR 
SED. For reference, SPIRE confusion limits on the HerMES 
fields are 5.8, 6.3 and 6.8 mjy beam " 1 at 250, 350 and 500 
fim, respectively dNguven et al.ll2010l ). In Table|4]we list the 
measured IR flux densities for the CO detected sources and 
candidates, including the radio source ID 51130. 

For the 3 galaxies that were detected at both 250 and 
350 /im, we fitted a single-component gray-body dust model 
using the Herschel photometry. We assume optically thin 
emission (r„ < 1), with the flux density given by: 



S„ = 



Dl 



M d K„[B„(T d ) - B V (T BG )], 



(1) 



where B V {T) is the Planck function, and and corre- 
spond to the dust temperature and mass. Dl and Tbg are 
the luminosity distance and cosmic microwave background 
temperature at the source's redshift, respectively. The dust 
absorption coefficient is parametrized as k„ = k,q(v/vq)^, 
where we adopt an emis sivity index /3 = 1.5 dKovacs et alj 
l2006l;lMagdis et al.ll201lD and nn = 0.4 cm 2 gr _1 at 250 GHz 
i Kruegel fc Siebenmorgenll 19941 ). 

Figure [13] shows the best fitting models to the far-IR 
data points along with the likelihood distribution, illustrat- 
ing the uncertainties involved in temperature and mass. We 
find dust temperatures of 34 ± 4 K, 30 ± 2 K and 33 ± 4 K, 
and dust masses of (2.3 ± 0.4) x 10 s M e , (7.3 ± 0.8) x 10 s 
M and (2.8 ± 0.5) x 10 8 M Q for ID 51613, ID 51858 and 
ID 51130, respectively. These correspond to total IR lumi- 
nosities (integrated over the range 8-1000/mi) of 5.7 x 10 11 
Lq, 9.5 x 10 11 Lq and 5.5 x 10 11 Lq, respectively. We note 
that given the lack of detections at shorter IR wavelengths 
(25 — 70 /xm rest-frame), it is not possible to estimate the 
contribution of a hot dust component, although the contri- 
bution to the IR luminosity is likely negligible. 

In the two cases where only the 250 /mi emission was 
detected, we fit only the dust mass (and hence the IR lu- 
minosity) by fixing the dust temperature to 35 K, which 
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Figure 13. Top: Best fit single-component gray body dust models for objects detected in two bands with Herschel. Bottom: Normalized 
likelihood distribution for different combinations of dust temperatures and masses for j3 = 1.5, obtained from the comparison between 
the data and the gray-body models. Contours represent levels at 0.4, 0.6, 0.8 and 0.9. 



is the ty pical value found for star-forming galaxi es at high- 
redshift (|Kovacs et all 120061 ; iMagdis et alThoill ). We find 
dust masses of 1.7 x 10 8 M and 1.1 x 10 s M©, and IR lu- 
minosities of 4.9 x 10 11 L Q and 3.3 x 10 11 L for ID 51207 
and ID 51380, respectively. Varying the assumed dust tem- 
perature in a typical range 30 — 40 K, leads to dust masses 
in the range (1.0-3.4) x 10 s M @ and (0.7-2.3) x 10 8 M Q , 
and thus IR luminosities in the range (4.2 — 6.0) x 10 11 Lq 
and (2.8-4.0) x 10 11 Lq, for ID 51207 and ID 51380, respec- 
tively. This corresponds to ~ 30% uncertainty in luminos- 
ity caused by our assumption in dust temperature. Finally, 
we remark that given the lack of sufficient constraints on 
the far-IR SED for these galaxies, our luminosity estimates 
could be affected by different choices of /?, as well as dust 
temperatures, and thus should be taken with caution. 



3.7 Galaxy properties 

We use the full multi-wavelength dataset available in the 
COSMOS field, including 31 broad, intermediate and narrow 
bands from the UV/optical and the IRAC 3.6^im, 4.5^im, 
5.8/im and 8.0/^m bands, to measure the range of physical 
properties of the observed galaxies in the JVLA field. 

We computed SFRs and stellar masses usi ng the SED 
fitting routine HyperZ ( Bolzoncll a et al.l |2000| ) and fixing 
the spectroscopic redshift as given by the optical spectra 
or the CO line identification, under the assumption that the 
latter is correct for the two CO identified candidates (ID 
51207 and ID 51380). We made use of an updated version of 



Table 4. Flux densities measured with Spitzer and Herschel 



1D a,b 


S 2 4 c 


S250 d 


5 35 o d 


Ssoo d 




(MJy) 


(mjy) 


(mjy) 


(mjy) 


51130 


220 ± 30 


9.7 ± 1.7(2.8) 


9.6 ±2.1(3.6) 


< 8.8 


51613 


125 ± 20 


10.4 ± 1.5(2.7) 


6.9 ±3.7(4.7) 


< 10.3 


51858 


150 ± 15 


19.0 ±2.6(3.4) 


18.7 ±2.3(3.7) 


< 10.4 


51207 


110 ± 15 


7.7 ± 1.5(2.7) 


< 9.6 


< 7.7 


51380 


65 ± 17 


5.0 ±0.5(2.3) 


< 13.3 


< 9.7 



"COSMOS ID; b Herschel PACS measurements led to 3cr upper 
limits of 6 and 12 mjy at 100 and 160 /im, respectively; c Flux 
density measured with the Spitzer MIPS at 24/^m; d Flux 
densities measured with Herschel SPIRE at 250, 350 and 500 
lira and were com puted using the 24/^m position as prior 
llLutz et al.ll201lh . Quoted uncertainties correspond to the map 
flux error and total uncertainty (including instrumental errors), 
shown in parenthesis, at the given band. 



the COSMOS photometric catalog fsee lSalvato et al.llioill ); 
and applied the appropriate galactic extinction corrections, 
and "auto-offsets" for SE D fitting computation (|llbert et al.l 
l2009l ; ISalvato et aDl201ll ) . 

We utilized 6 synt hetic template sets from the 
iBruzual fc Charlotl l|2003h libraries, including a constant 
star formation history (SFH) and 4 exponentially declining 
ones, with typical decaying timescales of 0.1, .3, 0.6, 3.0 and 
10 Gy r. We used an extinction correction from lCalzetti et ail 
(2000), with an optical extinction Ay = — 3. Both the best 
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Table 3. Physical properties 



ID a 


z b 

•^spec 


SFH t c 
(Gyr) 


Age d 
(Gyr) 


My 6 
(mag) 


A v f 
(mag) 


A/ S tars 3 
(Xl0 10 Mq) 


SFRsed h 
(Mq yr" 1 ) 


SFRjr i 
(Mq yr" 1 ) 


SFE i 
(Mq yr" 1 ) l 


50480 
51130 
51613 
51858 


1.523 
1.519 
1.516 
1.556 


0.1 

0.3 
0.6 
3.0 


0.25 
0.50 
0.36 
0.36 


-21.1 
-24.0 
-22.4 
-22.7 


0.2 

1.3 
1.9 
2.0 


o c + 0.2 
O.O_ 2 

24+ 4.0 

4 5+1-8 

6 0+ 16 
D - u -i.o 


4 
25 
130 
218 


94 
97 
162 


> 4 

> 80 
23 
75 


51207 
51380 


1.530 
1.551 


0.6 
0.6 


1.00 
0.25 


-21.9 
-21.4 


2.5 
2.6 


n +5.7 
11 -1.3 

2.8+H 


66 
120 


83 
56 


47 
29 



a COSMOS ID; b CO-based spectroscopic redshift or optical, if CO undetected ; c Timescale of exponentially declining star-formation 
history (SFH) best-fit template. d Age of the best-fit template SED; e Optical absolute magnitude (AB); f Best-fit extinction 
correction; 9 Best-fit stellar mass. Uncertainties correspond to the 90% confidence interval; h Star formation rate (SFR) derived from 
SED fitting; 1 Star formation rate derived from IR luminosity, assuming SFR/(A/q yr~ 1 )= 1.7 X 10" 10 Ltr,/-Lq; j Star formation 
efficiency, SFE= Lir/L^q in units of Iq = Lq (K km s _1 pc 2 ) -1 . 



fit template and the fixed redshift give the absolute magni- 
tude, which combined with a prescription for the mass-to- 
light ratio and SFR from the UV/optical continuum for each 
template, can be used to derive the stellar masses and SFR 
for each galaxy. Due to degeneracies in the extinction cor- 
rection and best template fit, and the use of templates that 
do not include emission lines, these estimates lead to typi- 
cal uncertainties of 30 — 50%, similar to what was reported 
in ot her similar studies (D addi et al. I l2010al ; iTacconi et alj 
l2010t ). Table [3] lists the obtained parameters for the best 
fitting template SEDs. For the CO detected and CO candi- 
date galaxies, we also show the best fitting template at the 
derived CO redshift (Fig. [14]). 

From the good agreement between the spectroscopic 
and the photometric rcdshift s for the gala xies in our field, 
we use the results from lllbert et all §009) as a cross-check 
for the parameters used in our fitting routine (e.g. template 
libraries, extinction). Such measurements provide an inde- 
pendent estimate of the stellar masses and SFRs at the best 
fit photometric redshift. We find that while both SFR esti- 
mates are in agreement with a typical uncertainty of 30%, 
the stellar mass estimates are in agreement to within 15%. 
The larger difference between both SFR estimates, com- 
pared to the difference in stellar masses, reflect the uncer- 
tainties in the extinction correction applied. While estimates 
for the stellar masses largely rely on the rest-frame optical 
and near-IR emission, which is less affected by extinction, 
estimates for the SFR are based on the rest-frame UV for 
which a proper extinction correction factor plays a key role. 

An independent and more direct measure for the SFR 
is given by the far-IR emission, and we use it as the ac tual 
SFR in the following of the text. From lKennicuttJ (|l99St ), we 
have SFR/lMQyr' 1 ] = 1.7 x 10~ 10 L m /L Q . The IR-derived 
and SED fitting derived estimates are in agreement to within 
25%, and have a dispersion of only ~ 10% from the average 
SFR (see Table EJ). The only exception is ID 51130, where 
the SED fitting routine implies SFR= 25 Mq yr -1 , while the 
IR luminosity indicates SFR= 93.5 Mq yr" 1 . This difference 
could be explained by the fitted template SED that does not 
represent the real nature of this source, despite it provides a 
good fit to the photometric data points (low \ 2 )- This large 
difference may also reflect the uncertainties in the derived 



IR luminosity given that the IR data is close to the SPIRE 
confusion limit. Table [3] lists the computed SFR values. 



4 DISCUSSION 

4.1 Star- formation rates and efficiencies 

Recent studies suggest that, compared to starbursting 
ULIRGs and SMGs, spiral and disk galaxies hav e substantial 
differe nc es in SFEs and dyn amical timescales (|Daddi et al.l 
l2010al fhl iGenzel et al.ll2010l ). This has been interpreted as 
evidence for different modes of star-formation. Disks appear 
to evolve secularly, with long dynamical and gas consump- 
tion timescales, while starbursts seem to evolve rapidly, pos- 
sibly triggered by merger-induced activity. However, it has 
been argued that such fundamental differences in the SFEs 
are, at least partially, produced by observational reasons. 
Aside from the large uncertainties in CO and IR measure- 
ments in ultra-luminous IR galaxies (ULIRGs) and SMGs, 
the large majority of the CO observations have been done in 
high- J CO transitions (J > 2). Extrapolation from high-J 
CO luminosities to CO 1 — 0, a direct tracer of molecu- 
lar gas mass, is not straightforward in a case by case basis, 
as CO 1 — measurements in a few SMGs found ~ 2x 
more mol ecular gas than pred i cted by the CO 3 — 2 emis- 
sion line l|Bothwell et alj|201fj llvison et ai1l201ll ; iRiechersI 
l201ll ). Furthermore, such differences in SFE are increased 
when using different CO luminosity to gas mass conversion 
factors for starburst and disks. Until accurate measurements 
of the CO 1 — emission line over a statistical number of 
galaxies (both disks and SMGs), it will be difficult to solve 
this issue. 

Although it is out of the scope of this paper to 
study these important differences, it is instructive to com- 
pare our CO 1 — measurements with other galaxy pop- 
ulations. Figure [15] compares the CO and IR luminosi- 
ties of different galaxy populations including our detec- 
tions in the cluster field. We do not convert CO luminosi- 
ties to molecular gas masses given the large potential un- 
certainties related to aco conversion factor across differ- 
ent galaxy populations. For the local Univer se, we com- 
piled measurements of typical spiral galaxies (|Lerov et al.l 
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Figure 14. Optical to mid-IR spectral energy distribution (SED) toward both CO detected galaxies (ID 51613 and ID 51858) and CO 
candidates (ID 51207 and ID 51380). Data points are drawn from th e COSMOS field multi-wavelength photometry. Solid lines show the 
best fitting templates from the libraries of Bruzual & Chariot (2003). 



l200Sl ) and local ULIRGs l|Solomon et al.lll997T i. Recent CO 
measurements of massive star-forming galaxies a t z ~ 



1.5 

- 2 



Daddi et al.l 



Tacconi et al 



2010al: 



l201Ch 



0.4 riGeach et al.l l201lt ). 
lAravena et al.l 20101 ') and z 
are also shown, as well as CO measurements of four lensed 
star-forming galaxies at z — 2 — 3 th at have intrinsic prop- 
erties similar to z ~ 1 — 2 disks dRiechers et al.l [2010b; 
iFraver et al.l l201ll ) . Additionally, we have included obser- 
vations of SM Gs and IR luminous QSO host galaxies from 
the l i terature (|Greve et al.l 120051: ISolomon fc Vanden BoutJ 



2005l;llono et al.ll200fj;IWeifi et al.ll2q07l; lAravena et al.l200St 
Coppin et al.l I2OO8I; IFraver et all l2008j: ISchinnerer et al 
20081 ; llvison et al.l l201ll; ICoppin et all I20T0I; ICarilli et al 



20ld. l201ll; iRiechers et alj l200d. l2010al ; iRiechersI l201ll : 
Riechers et al.ll2011al lbl; lFu et al.l l2012f r l 



1 In the case of the lGreve et al. I j2005h SMG objects, we recom- 
puted the IR luminosities by considering a optically thin, single 
component gray-body dust model of the form S„ oc BuiT^M^, 
with a dust temperature of 37 K and /3 = 1.5, typical of star- 



The SFE can be defined as SFE= Lib./L' co in units of 
Lq (K km s _1 pc 2 ) -1 . For galaxies detected with Herschel, 
we directly use their estimate of the IR luminosity. For those 
galaxies with no Herschel detection, we adopt the SED fit- 
ting derived SFR and convert them int o IR luminosities fol- 
lowing the relation of iKennicuttl |l998). The obtained SFE 
values are shown in Table [3] For the CO detected galaxies, 
the SFE ranges between 20 to 75 L Q (K kms" 1 pc 2 ) -1 . 

Figure [TBI shows that the SFE from the galaxies in our 
sample are overall consistent with the general population 
of star-forming galaxies at high-redshif ts, which range be- 
tween 30 to 180 L Q (K km s" 1 pc 2 )" 1 (|Daddi et alj|2010al : 



forming SMGs l lKovacs et alj|2006h . constrained by the source's 
redshift and 350 GHz flux densi ties. This implies I R luminosi- 
ties ~ 2x lower than reported bv lGreve etal] <200Sl) . For detec- 
tions made only in a few high- J CO transitions (J > 1), with no 
CO excitation models available, we ado pt line temperature ra- 
tios: R31 ~ R41 ~ 0.5 and R21 ~ 0.8 (Ivison et al. 2011; R iechers! 
120111) . 
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Figure 15. CO J = 1 — and IR luminosities of the CO-detected galaxies in our clust er field, shown by black filled circles, compared 
with other gala xy samples. Black c rosses show measurements for local spiral galaxies iLerov et al.ll200St) and green asterisks represent 
local ULIRGs {Solomon ct al.|[l997t) . Grey triangles and squares show high -redshift SMGs and QSO hos t galaxies, respectiv ely (see text). 
Red triangles, squares and circles show m assive disk galaxies at z ~ 0.4 iGeach et al 1 l201lf) . 2 ~ 1.5 (Da ddi et a"l]|2010al') , and typical 
massive star forming galaxies at z ~ 1 — 2 llTacconi e t al. 2010), respectively. Red 5- pointed stars show four lensed massive st ar-forming 
galaxies at z = 2 — 3, which have intrinsic properties similar to z ~ 1 — 2 disks i Riechers et alj|2010bl: i Fraver et al.ll2011f) . In cases 
were only J > 1 CO lines are available, we corrected to the CO J = 1 — line (see text). Following iDaddi et alj d2010al 1 . we show a 
representative fit to local spirals and disk galaxies at high-redshift (grey dashed line), logio Lir = 1.12x logioLco + 0.5, for guidance. 
The dotted line shows the same line, with a factor +0.5 added. 



lAravena et al.|[2010| ; iTacconi et al]|2010l ). however, the sig- 
nificant uncertainties in our measurements make it difficult 
to draw further conclusions. On the other hand, the agree- 
ment between the SFE of ID 51207 and ID 51380 with typi- 
cal high-redshift galaxies, gives fundamental support to the 
identification of those objects as CO line emitters. 

As mentioned before, the optical source ID 51861 has an 
optical spectroscopic redshift of 1.240 although it is wrongly 
identified with a line peak. In the absence of a spectroscopic 
redshift, it would be identified with CO 1 — at z = 1.497. 
This source is not detected at 24/im and hence, by definition, 
not detected with Herschel since we are using the 24 /im 
prior Herschel catalog. It is interesting to note that if we 
had identified this galaxy as a CO emitter at the respective 
redshift, the SFR obtained from SED fitting of 10 M Q yr" 1 
would have implied a SFE of ~ 4 L (K km s -1 pc 2 ) -1 . This 
would have been lower than the ones measured in local spiral 
galaxies, wh ich are typically in the range 10—100 Lq (K km 
s" 1 pc 2 ) -1 (|Lerov et alj|2008h . yet still consistent with that 
found in local ellip tical galaxies, in the range ~ — 50 Lq 
(K kms" 1 pc 2 )" 1 (jShapiro et al.ll20irJ ; [Crocker et al.ll201lD . 
If we did not have an optical spectroscopic determination of 



the redshift, such mis-identification would have still been 
suspicious taking into account such low SFE and the overall 
multi- wavelength properties of the host galaxy. 

4.2 Comparison to the predicted CO luminosity 
function at z — 2 

One of the major motivations to perform blind deep CO sur- 
veys of the sky is to measure possible evolutionary trends 
in the amount of molecular gas with redshift and CO lu- 
minosity (or gas mass). An important tool to do this is 
to estimate the CO luminosity function at different red- 
shifts. Combining results from the Millennium simulation 
( Spr ingel et al. I l2005l ; iDe Lucia fc Blaizotl [20071 ) and empir- 
ical, detailed models of abunda nces of molecular gas, an d 
their relation to CO intensities, lObreschkow et all (2009), 
made predictions for the luminosity function of the CO 
emission lines as a function of redshift and CO rotational 
transitions. In this section, we measure the space density 
of CO detected galaxies within the physical volume covered 
by our observations, and co mpare this with the value s ex- 
pected from the predictions l|Obreschkow et aT]|2009l ) and 
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Figure 16. Predicted cosmic evolution of the space density of 
CO 1 — emitting galaxies as a function of CO 1 — luminos- 
ity (CO luminosity function) at redshift s z = (gray solid line) , 
2 = 1 (black dashed line) and 2 = 2 l lObreschkow et al J 12009 . 
black solid line; ). The gray open triangles show th e space den- 
sity o f CO emitters measured in the local Universe ijKeres et alj 
2003). The black filled circle shows the space density derived from 
our 4 CO detections in the overdense field targeted with JVLA. 
The gray filled circle shows the measurement but decreased by 
7 to correct for the overdensity. The black and gray open dia- 
monds represent the space density obtained considering only the 
2 spectroscopically confirmed galaxies uncorrected and corrected 
by the 7x factor, respectively. The black open square shows the 
estimate of space density of CO emitters derived from the 6 BzK 
star- forming galaxies at z ~ 1. 5 jDaddi et al.ll2010al 'l. 



with previ ous observations of massive star-forming galaxies 
atz~ 1.5 (|Daddi et al.ll2010al ). 

To estimate the space density implied by our 4 CO 
detections, we start by computing the CO luminosity in 
units of Jy km s" 1 Mpc 2 as L v = (1 + z)~ 1 4nDtS v 
(adopting the terminolo gy presented in the appendix A of 
lObreschkow et all [20091 ) , where z is the redshift, Dl is 
the luminosity distance in Mpc, and S v is the velocity- 
integrated CO intensity in units of Jy km s~ . From our 
4 detections, we estimate a weighted average CO luminos- 
ity of (6.3 ± 1.6) x 10 7 Jy km s" 1 Mpc 2 . Next, we deter- 
mine the cosmological volume covered by our observations, 
considering the area comprised by our JVLA field of view 
(» 60" x 60") and the redshift coverage, z = 1.478 - 1.591. 
This leads to a volume of AV = 356 Mpc 3 . With this, we 
compute the space density of CO sources in units of Mpc -3 
dex -1 as cj>o = N/(AVlog 10 Lco), where N = 4 is the num- 
ber of CO detected galaxies (Fig. 1161) . 

We also compute the space density of CO emitters im- 
plied by the sample of 6 massive star-forming BzK galax- 
ies _tha^_werereUably detected in CO 2 — 1 line emission 
bv lDaddi et al. | (|2010al ). These are typical massive galaxies 
with an average IR luminosity of ~ 2 x 10 12 Lp. For this, 
we fo llow the procedure used by lObreschkow fe Rawlingsl 
(2009), but we adopt more recent developments to our anal- 
ysis, including the use of 6 BzK galaxies available instead 



of only 2 from iDaddi et al.l |2008l ). For each object, we 
compute its CO 2 — 1 luminosity directly fro m its redshift 
and me asured CO line intensity presented by IDaddi et al.l 
l|2010al ). We convert the CO 2 — 1 luminosities (in Jy km 
s _1 Mpc 2 ) into CO 1 — luminosities by assuming that 
both transitions are typically in local thermodynamic equi- 
librium (LTE) i n BzK galaxies , and t herefore scaling approx- 
imately as v 2 (jAravena et al.|[2010l ). Thereby, we compute 
the weigh ted average CO 1 — luminosi ty for the sample. 
Following lObreschkow fc Rawlingsl (|2009l ) , we approximate 
the space density of CO sources as the space density of IR 
sources at Lib. ~ 2 x 10 12 Lq, which corresponds to the 
average IR luminosity for all 6 objects. For typical galax- 
ies at this IR luminosity level a nd in the redshift range 
2 = 1.2 — 2.0, Rodighic ro et al.l (|2010l ) finds a space den- 
sity $ir = 3.2to.5 x 10" 4 Mpc -3 per unit \og w L m , which 
we scale to logiol/co to find the space density of CO sources 
(Fig. 1161) . We selected only this sample of massive star- 
formi ng galaxies at high r edshift, with respect to the sample 
from lTacconi et al] |20ld ). because (1) all of them were de- 
tected in CO 2 — 1 with high significance; (2) they have 
similar redshifts compared to our CO sources; and (3) three 
of these sources correspond to the only 'normal' objects so 
far detected in CO 1 — emission, making a determina- 
tion of line ratios between the 2 — 1 and 1 — transitions 
more direct and straightforward. We remark, however, that 
this method to estimate the space density of CO emitters 
is limited by two important assumptions: the main one that 
the space density of CO emitters is directly traced by the 
space density of IR sources at a given IR luminosity does 
not take into account that there appears to be different star 
formation laws for bright starbursts and the general popula- 
tion of galaxies (|Daddi et al.ll2010bl . see previous section;). 
Without knowledge of the nature of individual sources in 
the IR sample from where we adopt the space density, it 
is difficult to estimate the uncertainties introduced by this 
assumption. Furthermore, the 6 BzK galaxies were selected 
to be massive, K-band selected, galaxies with significant 1.4 
GHz radio detections, and thus not strictly IR-selected as 
it is assumed above although all 6 sources share similar IR 
luminosities. The actual density of such BzK sources is thus 
likely lower than the space density of IR sources in the field. 
Hence, this procedure should be taken with caution and only 
for comparison purposes. 

Figure \W\ shows the space density of CO sources as a 
function of the average CO 1 — luminosity obtained for our 
covered field as well a s the space d e nsity o f CO sources from 
the 6 BzK galaxies of Daddi et al. (2010a), compared to the 
predictions from lObreschkow et al.l (|2009l ) for z = 0, z = 1 
and z = 2 and observations of lo cal galaxies by Keres et al.l 
(2003). Similar to the findings of lObreschkow et al. I (|2009l ). 
the space density of CO sources derived from the BzK galax- 
ies is only marginally consistent with the one predicted by 
the simulations. Furthermore, the space density obtained 
for the CO sources in our field appears to be in excess from 
that predicted in the simulations and from that obtained 
for the BzK sample. This result comes naturally from the 
fact that our observations likely target an overdensity of 
galaxies at z ~ 1.5. If we correct the ~ 7x increase factor 
due to the overdense field, we obtain a value in agreement, 
within the uncertainties, with what is measured from the 
BzK galaxies and marginally consistent with expectations 
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from simulations. Similarly, if we take into account only the 
two spectroscopically confirmed CO sources in our field, the 
space density and average CO luminosity tend to agree bet- 
ter with the estimate obtained from the BzK galaxies. The 
apparent offset between measurements and the predictions 
from simulations could be explained by several reasons or 
a combination of them, for instance: (a) lack of a numer- 
ous sample of CO sources; (b) although the Obreschkow et 
al. model assumes a galactic CO luminosity to gas mass 
conversion factor, the use of different values could produce 
a better agreement between predictions and observations; 
and (c) source-to-source variations in the line ratios, e.g. 
between 2—1 and 1 — 0. These results exemplify the need 
to perform measurements of CO emission, particularly CO 
1 — 0, in a significant sample of galaxies (ideally in a contigu- 
ous field) selected based on their 1R luminosities, in order to 
constrain the evolution of the CO luminosity function and 
thereby models of galaxy formation and evolution. 



4.3 Future searches for blind CO line emission 

The determination of redshifts for galaxies based on 
(sub) millimeter and radio wavelength observations have 
been focused on the 1R brightest sources, favored by the 
expanded bandwidth of such facilities. Based on the re- 
lationship between IR and CO luminosities, an IR bright 
galaxy, such as an SMC, is likely to also have bright 
CO emission. This facilitates a high significance "blind" 
CO emission line detection along the line of sight, lead- 
ing to an accurate redshift measurement. In such stud- 
ies, the position of their targets are known a priori, usu- 
ally derived from large blank- field IR/submm surveys, per- 
mitting to spatially guide the CO observations, bein, 
thus not really 
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Combes et al.ll2012 ). However, serendipitous CO line detec- 



tions, where the locations and redshifts of the g as-rich galax- 
ies we re not anticipated hav e been reported ([Daddi et alj 
l2009bl lal; iLestrade et aLll201Gft . 

Blind spectroscopic imaging of the sky with large band- 
widths have a huge potential for discovery of gas-rich, op- 
tically and IR faint galaxies, however they require deep in- 
tegrations to ensure significant emission line detections, at 
the ~ lOer level, and unequivocally identify the galaxies and 
their redshifts. In the absence of such high signal-to-noise 
detections, observations need to be guided by well-known 
optical, IR or radio positions and by accurate photometric 
redshifts. As shown in this work, focused line observations 
of clustered fields at high-redshift have the advantage to in- 
crease the likelihood of finding gas emission using a limited 
bandwidth. This approach reduces the required observing 
time compared to purely blank fields, where besides the need 
to cover a substantial area of the sky, it is necessary either 
to cover a reasonable range in frequency (e.g. a full 3 mm 
band for ALMA, a full K or Q band with JVLA) in order 
to contain emission lines at different redshifts, or to perform 
"spot-checks" at a given redshift. 



5 SUMMARY AND CONCLUDING REMARKS 

In this work, we have presented deep CO 1 — line obser- 
vations of galaxies located in a galaxy cluster candidate at 
z — 1.5 using the JVLA. The candidate cluster was iden- 
tified using the red-sequence technique and is associated to 
an overdensity of ~ x7 of star-forming galaxies. We use the 
spatial clustering and expanded bandwidth of the JVLA to 
simultaneously observe the CO emission from four galax- 
ies with available optical spectroscopic redshifts in a single 
pointing and frequency setting. 

We detect, at the ~ 3cr and ~ 4<r level, the CO 1 — 
emission line from two of the galaxies with available spec- 
troscopic redshifts: ID 51613 and ID 51858 at z — 1.516 and 
z = 1.556, respectively. We find L' co = (3.9 ± 1.1) x 10 10 K 
km s" 1 pc 2 and L' co = (1.3 ±0.4) x 10 10 K km s _1 pc 2 , re- 
spectively. Both galaxies are also detected with Spttzer at 24 
/im, and with Herschel at 250 /im and 350 /im. While model- 
ing of their dust properties indicates total IR luminosities of 
5.7 x 10 11 Lq and 9.5 x 10 11 Lq, respectively, characteriza- 
tion of their multi-wavelength SEDs suggests both are young 
star-forming galaxies that have formed a significant fraction 
of their stellar content, with stellar masses of 4.5 x 10 10 Mq 
and 6.0 x 10 10 Mq, respectively. 

We performed a blind search for significant (> 4a) emis- 
sion line peaks in the JVLA data cube around the posi- 
tion of optical sources that were selected to have photo- 
metric redshift in the range z — 1.2 — 1.8 and limited to 
Ks < 23.5. This selection criteria ensures an accurate pho- 
tometric redshift, while selecting galaxies with stellar masses 
> 2 x 10 9 Mq . We find that 2 of these selected optical galax- 
ies are associated with significant emission line peaks (> 4a), 
which are thus identified as CO 1 — line emission candi- 
dates. Both galaxies, ID 51207 and ID 51380, are found to 
have photometric redshifts in the range 1.4 — 1.6. Based on 
Monte Carlo simulations of the distribution of such emission 
line peaks in the field, we find that any optical galaxy that 
follows our selection criteria and is located within 1" from 
a candidate emission line peak is statistically unlikely to be 
associated to such peak by chance. In the case of ID 51207, 
such probability is 3.7% while for ID 51380 we find 10%. 
Using in parallel a Bayesian inference approach to select 
emission line peaks in the data cube, we recover both CO- 
emitting candidate sources. Here, the line identified with ID 
51380 is found to have a probability of 70% of being a CO 
emission line, however the line identified with ID 51207 is 
found to have a rather low probability of 13%. 

We find CO based spectroscopic redshifts of 1.530 and 
1.551 with L' co = (1-03 ± 0.25) x 10 10 K km s _1 pc 2 and 
L' co = (1-12 ± 0.26) x 10 10 K km s' 1 pc 2 for ID 51207 
and ID 51380, respectively. These galaxies, ID 51207 and 
ID 51380, are detected at 24 /im and 250 fim with Spitzer 
and Herschel. Such measurements lead to IR luminosities of 
4.9 x 10 11 Lq and 3.3 x 10 11 Lq for ID 51207 and ID 51380, 
respectively. In the former case, its optical/IR SED suggest a 
~1 Gyr old, very massive star forming galaxy, with a stellar 
mass of 1.1 x 10 11 Mq, while for ID 51380 the SED indicates 
that it corresponds to a young star-forming galaxy, with a 
stellar mass of 2.8 x 10 10 Mq. 

In one case, we associated a line peak candidate with 
an optical source that had a previous optical spectroscopic 
redshift. This secure optical redshift of z — 1.240 is incom- 
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patible with the one that would be derived from the identifi- 
cation of the line peak with the CO f — emission line. This 
optical source fits well within our source selection criteria, 
K < 23.5 (AB) and z p hot = 1.2—1.8, however if we had iden- 
tified it with the CO line emission, the non-detection in the 
1R bands and its implied SFR from the SED fitting method 
would have implied a very low SFE. We argue that in the 
absence of high-significance line detections, tight restrictions 
should be imposed on sources to be identified with CO line 
emission, such as detection in the IR bands or a lower limit 
in SFRs. 

We measured the space density of CO galaxies com- 
pared to the space density of CO emitters estimated from 
the 6 BzK galaxies significantly detected in CO emission 
by Daddi et al., and compared to predictions from semi- 
analytic simulations. Overall, we find that all observations 
are only roughly consistent with the simulations, despite the 
low number of detections of typical star-forming galaxies at 
high-redshift. Clearly, observations of molecular gas from a 
statistically significant sample of these galaxies is necessary 
to measure the evolution of CO luminosity function with 
redshift, and thus constrain models of galaxy formation and 
evolution. We conclude by discussing the advantages of per- 
forming deep CO observations of star-forming galaxies in 
clustered fields compared to blank-fields in the sky. 
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